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Abstract 

We have succeeded in applying Ziglin's test on non-integrability to non- 
homogeneous nonlinear lattices (Fermi-Pasta-Ulam lattices). By explicit cal- 
culations of the eigenvalues of the monodromy matrices concerning the nor- 
mal variational equations of Lame type and with the use of the non-resonance 
hypothesis about the eigenvalues, we obtained a theorem proving the nonexis- 
tence of additional analytic conserved quantities other than the Hamiltonian 
itself for FPU lattices in the low energy limit. Furthermore, after introducing 
a concept of degree of non-integrability, we have investigated the classifica- 
tion of non-homogeneous nonlinear lattices using a transformation 3? from 
a non-homogeneous nonlinear lattice to another non-homogeneous nonlinear 
lattice, which preserves their degree of non-integrability. 
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1 Introduction 



Recently, we know two major sources of detecting non-integrability of Hamiltonian sys- 
tems. One is Poincare's theorem|| of nineteenth century and the other is Ziglin's 
analysis |19| relatively recently found in the connection with the classical singularity 



analysis|9[] by Fuchs, Kowalevskaya and Painleve. What matters here is the following 
thing. The two different methods have significantly different predictive capacities to tell 
non-integrability. While Poincare's theorem is asserting the non-existence of additional 
conserved quantities such as 

$o({U) + fJ.$i({qu,Pu}) + At 2 $ 2 ({g„M) + • • • , (1-1) 
for a Hamiltonian system perturbed from an integrable systems 

H = H ({q u }) + nHriiq^Pv}) + fi 2 H 2 {{q v ,p v }) + ■■■, (1.2) 

where \l is a coupling constant whose explicit value is unknown, Ziglin's type theorem 
can prove the non-existence of additional analytic integral for an explicit Hamiltonian. 
To illustrate the differences between them more concretely, we can give the following 
examples. The restricted three-body problem was proved to be non-integrable by Whit- 
taker in 1904[Q based on Poincare's theorem, while Ziglin's theorem have never proved 
its non-integrability because of the problem that a Hamiltonian of the restricted three- 
body problem cannot be defined without the infinitesimal operation (m — > 0) and we 
don't have any explicit Hamiltonian which is necessary to perform solid Ziglin's anal- 
ysis. On the other hand, Poincare's theorem assumes the condition on the Hessian of 
Hamiltonian, i.e., the Hessian of Hamiltonian should not vanish 

V 0. (1.3) 



dq v dq x 



By this assumption, Poincare's theorem cannot not tell anything about the non- 
integrability of nonlinear lattices consisting of harmonic terms Hq = | Z)£=i Qi P ms 
unharmonic terms such as the nonlinear lattices of Fermi-Pasta-Ulam type (FPU lat- 
tices). Therefore, in the present paper, we examine the possibility of non-integrability 
proof of FPU lattices via Ziglin's analysis. It should be also mentioned that Ziglin's 
analysis have some technical difficulties, as is illustrated later, and no non-integrability 
proof has been so far studied either for realistic Hamiltonian systems such as the FPU 
lattices with non-homogeneous potential functions, except few rare instances such as the 
Henon-Heiles systems with two degrees of freedom^, ^(| |6| or for Hamiltonian systems 



with many degrees of freedom besides the case of Hamiltonian systems with global and 



symmetric coupling 13 



The purpose of the present paper is to show that this Ziglin's analysis can be per- 
formed for checking the non-integrability of nonlinear lattices including the FPU lattices 
and to show that these nonlinear lattices have no other analytic conserved quantities. 

In Section 2, we show that singularity analysis can be done for a class of FPU lattices 
by giving the explicit eigenvalues of the monodromy matrices about a special solution in 
terms of elliptic functions. 
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In Section 3, we prove a non-integrability theorem for FPU lattices based on the 
singularity analysis in Section 2, which shows that each FPU lattice has no additional 
analytic conserved quantities in the low energy limit. In Section 4 we discuss the classi- 
fication of non-homogeneous nonlinear lattices via the degree of non-integrability. 

In Section 5, we give a brief summary and open questions about the present analysis. 

2 Singularity analysis of nonlinear lattices 

Consider the following one- dimensional lattice: 

i n n+1 

z i=i i=i 

If this Hamiltonian admits the reflective invariance, namely, if it is invariant under the 
involutive symplectic diffeomorphism, 

J r H = H, (2.2) 

where 

J r : C 2n -> C 2n , J r : (gi, • ■ ■ ,q n ,Pi, • • • ,p n ) (-gi, • • • , ~q n , ~Pi, • -Pn), (2.3) 

the potential function v(X) must satisfy the condition given by 

v(X)=v(-X). (2.4) 

In the case of harmonic lattices, we have v(X) = \\iiX 2 . Thus more general nonlinear 
lattices with the reflection symmetry can be defined by the following polynomial function 
of even degree in X: 

V ( X ) = ^X 2 + ^X 4 + • • • + ^±X 2m . (2.5) 
v ' 2 4 2m y ' 

In this paper, however, we assume that this potential function is a polynomial function 
of degree 4, namely, we consider the Fermi-Pasta-Ulam(FPU) lattices^ 

in n+1 . n+1 

Hfpu = - 9 Y,P* + ?E(*-i - + f Efe-i " g,) 4 - (2-6) 

Z 1 = 1 Z 1=1 4 8=1 

Imposing the fixed boundary condition as 

go = gn+i = 0, n = odd, (2.7) 
the equations of motion are given by 

gi = -M 2 gi - g2) - Mg? + (gi - ga) 3 ), 

g2 = -/i2(2g 2 - gi - qa) - ^((ci2 - gi) 3 + (92 - g3) 3 ), 

(2.8) 

q n = -A*2(2g„ - g„-i) - yu 4 ((gn - g«-i) 3 + g 3 ) 
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Define M' = {(q,p) e C 2n \q x ^ or p 1 ^ 0}, M = M'/J and the canonical 
projection 7r : M' — > M. An additional complex analytic integral in C 2n , if it exists, also 
induces the corresponding analytic integral on M. Therefore it is sufficient to check the 
non-integrability of the reduced equations of FPU lattices on M in order to confirm the 
non-integrability of the original equations of FPU lattices. 

Consider the following one-parameter family of the solution for (|2.8| ) 



T(e,t) : q 3 = (( ^ l \ -\)^C<t>{t) (2.9) 
with the initial condition 

0(0) = 1, 0(0) = 0, H FPU (0) = e. (2.10) 
This is a complex solution in the complex time plane as 

71— 1 

qi = C<f>{t), q 2 = 0, q 3 = -C(f>{t), ■ ■ • , g n _! = 0, q n = {-!) — C4>(t), (2.11) 

which is depicted in Fig. 1. Thus, the underlying equation of <fi(t) is reduced to the 
following Hamiltonian system with one degree of freedom: 

4> + 2/i 2 + 2/i 4 C 2 3 = 0, (2.12) 

whose Hamiltonian is given by 

H{il>,4>) = \(<P) 2 + ^ 2 + ^yV = Const. (2.13) 
Then the original total energy e is written as follows: 

e = H = H(d>, ft^C 2 = ^C\^ 2 + \^C 2 ) (2.14) 

by the initial condition ( 2.1C ). This equation ( |2.14| ) enable us to choose C as the following 
complex quantity: 



C 

Therefore, H(<b,(b) becomes 



\ 



A + ^-jh (215) 



H{q,p) = -p +/i 2 g +~ q, (2.16) 

where q = <f),p = <fi. By combining ( |2.13| ) with ( |2.14| ), we obtain the differential equation 
of (p(t) as 

^(0) 2 = 7 2 (l-0 2 ) + y(l-0 4 ) (2.17) 

with the system parameters 

72 = // 2 , 74 = fi 4 C 2 . (2.18) 
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Thus, the phase curve T(t, e)(|2.9|) is given by the elliptic integral 



and the solution of this differential equation (|2.f 7|) with the condition 

72m=4 7^ 

is given explicitly by the formula 

4>(t) = cn(k; at), 

where 

a = J2j 2 + 274, k ~- 



7i 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



/ 272 + 274' 

cn(k; at) is the Jacobi cn elliptic function, and k is the modulus of the elliptic integral, 
because of the following equalities 



#) 2 = {-asn(k;at)dn{k;at)} 2 

= a 2 {(-2k 2 + 1)(1 - cn 2 (k; at)) + k 2 (l - cn\k; at))} 
= 2 72 (l-0 2 )+ 7 4(l-0 4 ) 



(2.23) 



in Eq. (|2.17 ). We remark here that the relation 

72 + 74 = A*2 + C 2 fi 4 = 



V2 + 



4f 



n + 1 



Hi > 0, 



(2.24) 



holds for /14 > 0, [12 > 0, e > 0. This means that the modulus of the elliptic function k 
satisfies the following relation: 



< k 



1 

V2 



\ 



1 1 4e/j 4 
"T (n+1)^ 



< 



V2 



(2.25) 



and C is a real quantity. On the contrary, if the following relations 

/j, 2 > 0, /i 4 > 0, e < 



(2.26) 



are satisfied, the modulus of the elliptic function k becomes pure imaginary, i.e., k 2 < 0. 
C becomes also pure imaginary as in Fig.l. We note here that even in this unphysical 
case, there is no barrier to perform the following singularity analysis. However, in the 
following we shall consider the physical case with the condition as \i 2 > 0, /X4 > 0, e > 0. 
According to the reduction tc : M' — > M = M'/J by the reflective symmetry J, the 
solution ( |2.11| ) of the induced FPU lattices on M for /i 4 > 0, fi 2 > is a single-valued, 
meromorphic and doubly periodic function in t, having the two fundamental periods in 
the complex time plane as follows: 



2K(k) 2K(k) + 2%K'{k) 

Ti(e,/i) = — : — , T 2 (e,/x)- 



a 



a 



(2.27) 



Here K(k) and K'(k) are the complete elliptic integrals of the first kind: 



K(k) 



dv 



v 2 ){l - k 2 v 2 ) 



K'(k) 



dv 



v 2 )(l 



k 2 )v 2 



(2.28) 



Note here that these solutions have a single pole at t = r, where r = 2K ^ +i — ^ (mod T 1; T 2 ) 
in the parallelogram of each period cell; the phase curves T(e,t) are punctured tori. If 
we consider the solution ( |2.11| ) not on the reduced phase space M = M' / J but on the 
original phase space M', the solution has the following fundamental periods as 



T , ( v mk) rpi i s 2K(k)+2iK'(k) 



a 



a 



(2.29) 



which are different from ( 2.27| ) and this solution on M' has two poles at t = r', 

where r' = 2K< ^ +% K ^ ; 4E "( fc ) _|_^ - ftr ( fc ) (mod T^, T 2 ) in the parallelogram of each periodic 

cell. ' 

Let us consider the variational equations along the phase curve T(e,t) in ( [2.9|) of 
the general nonlinear lattices with the reflection symmetry. The linearized variational 
equations are written as follows: 



vj = Zj 



En 
k=l 



d 2 V 



dq k dqj 



-(72 + 3 74 cn 2 (fc; or*))(2& - ^-i - ^+i) for 1 < j < n, 



(2.30) 



where £ = £n+i = Vo = Vn+i = and £ d = Sq^ iy = 5pj (1 < j < n). 

If we rewrite these linear variational equations ( |2.30|) in the form of the vector variational 

equation as 



d^ 
dt 2 ' 1 



-(72 + 37 4 cn 2 (/c; at)) 



2 -1 ■•• 
-1 2 -1 ... 
-1 2 -1 ... 



••• -1 2, 

and use the fact that the eigenvalues of the n x n symmetric matrix 



(2.31) 



G 



2-10 ••• 
-1 2 -1 ••• 
0-12-1 












-1 



(2.32) 



are obtained as {4sin 2 



■2(n+l). 



|1 < j < n} by a normal orthogonal transformation G 



OGO 1 , the variational equations ( |2.30| ) can be transformed into the decoupled form: 



-4sin^ 



"2(71+1 



0(72 + 3 74 cn 2 (fc; «*))£■(*) (1 < j < n), 



(2.33) 
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where £' = 0$,. Clearly, these equations are regarded as the vector form of Hill's 
equations^ 

^ + A{t)£' = 0, A(t + T) = A(t), (2.34) 
where T = T\, T 2 . For j = the relation 



&±i = \M-r(6 - 6 + & + • • • + (-l)^n) (2.35) 

2 y n + 1 
holds. The corresponding variational equation 

C'jtti = -2( 72 + 37 4 cn 2 (A;; aO)Ckl(*) (2-36) 

2 2 

has a time- dependent integral /(£, t) = 77; t) derived from the variational operation 
D applied to the Hamiltonian (|2.6|) : 



ri ;t) = DH=(r r ^ + $- ±)H = rj-p + ^-V q 

= C<j>( Vl - m + m + --- + (-l)^Vn) (2-37) 
+2{C l2 cn{k- at) + C^cn 3 (k; at))fa - £ 3 + £5 + • • • + (-I)tt^), 

because 



5S = *(6 -& + •■• + (-1)^6.) 



n-1 



(2.38) 



+20( 72 + 3 74 cn 2 (A;; at)) fa -& + ••• + (-l)^Cn) = 0. 

Eq. ( |2.36| ) is called the tangential variational equation. On the other hands, a (2n — 2)- 
dimensional normal variational equation(NVE) is defined by the equations ( |2.33| ) with 
the tangential variational equation fl2.36| ) removed as follows: 

V'j = -4sin 2 (^ Ty )( 7 2 + 3 lA cn\k; at))%, 
^ = r lj 

for 1 < j& s±l) < n. 

Let us consider the monodromy matrices g defined by the analytic continuation of 

the solution C'(t) = (£(*), ^(t), • " " > &±i(*),%±i (*), • • " > of the NVE (p|) 

2 2 

along the periodic orbits in the phase curves T(e,£) as follows: 

C'(Ti) = gxC(0), C'{T 2 ) = g 2 C'(0). (2.40) 

The periods in Eqs. (|2.40| ) are Ti,T 2 in Eqs. ( |2.27|) , respectively. These two funda- 
mental periods T\ and T 2 naturally form the parallelogram, whose associate monodromy 
matrices are given by gig29i l g 2 l {= 9*)- These monodromy matrices are naturally en- 
dowed with the symplectic structure and the pairing properties of the eigenvalues, namely 
{<7i, erf 1 , 02, (T 2 l , • • • , o~ n , cr" 1 } . In general, the explicit calculation of the eigenvalues of 
the monodromy matrices is an unsuccessful business except rare cases such as Hamilto- 
nian systems with homogeneous polynomial functions jL7|, Riemann's equation, and the 
Jordan-Pochhammer equations^Xl ^, E9. This is one of the unavoidable difficulties in 
performing Ziglin's analysis of general dynamical systems. In our case, however, we can 



7 



^| - (E lS n 2 (k; at) + E 2 )y = 0, (2.41) 



compute explicitly the eigenvalues of the commutator = gig 2 g± l g 2 1 by using the fact 
that the normal variational equation (|2.39|) happens to be in a class of Lame equations 


where E\ and E 2 are constants, and the eigenvalues a of the commutator g* = g\giQ\ l <h 
are known [15 to be determined by the indicial equation 

A 2 - A - (a 2 k 2 )E 1 = 0, a = exp(2niA) (2.42) 

with the singular point (pole) r located at the center of the parallelogram: 

= T 1 + T 2 = 2K{k)+iK\k) 

2 a 1 ' 1 

If we apply the indicial equation( |2.42| ) to the normal variational equations (2.39) of the 
FPU lattices, the exponents of the eigenvalues of g* are given by 

A 2 - A - 124^sin 2 ( , jn - ) = 0. (2.44) 



Noting 

74 



a 2 k 2 

from ( |2.22| ), we finally obtain the eigenvalues of the commutator as 



(2.45) 



1± Jl + 48sin 2 ^4TT 

/r> • V 2(n+l) 

Uj = exp{2m 



= -exp(±7riy25 - 24cos(^y)) (2.46) 

for 1 < j(^ — n -> because there is only one pole singularity inside the parallelogram 
formed by counterclockwise closed loop of the monodromy group g\g29\ X &2 ■ Even for 
more general nonlinear lattices given by potential functions ( |2.5| ) of degree 2m we can 
also compute the phase factors. Let us consider singular solutions <p(t) and £j(t) whose 
singularities are located at t = r in the complex time plane as follows: 

m=C\t-rf ) (3<0 

m = c"(t-rp, [ZA() 

where C is a particular constant, C" is an arbitrary constant. Thus, we can easily 
obtaine the formulas the relations for /3, C, uj as 



rra-l' 

(C') 2m - 2 = (2.48) 
, 2 -, J -2^)sin 2 (^ TI ) = 
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by the equations of motion. With the use of the fact that the phase factor is obtained 
by means of the analytic continuation along the closed loop around the singular point 
t v ' , we obtain the phase factors as 

exp(27uz/j) 

i± « A+8 m } 2m -y sin 2 
= exp{27rz— ^ ^i^)} (2.49) 

= -expiiTr^l + S^^sin^^)} 

where 2m is the degree of the potential polynomial of v(X) in ( |2.5|) . It is easy to confirm 
that in the case that 2m = 4, the formula (|2.49|) recover the eigenvalues Q2.46| ) of the 
monodromy matrices obtained from the Jacobi elliptic function. However, in the other 
case that 2m > 4, we cannot compute the eigenvalues of monodromy matrices only from 
the phase factors in general, because the normal variational equations do not belong to 
a class of Lame equations. To summarize, we can say that it is essential to perform 
the singularity analysis towards non-homogeneous nonlinear lattices using the explicit 
eigenvalues of the monodromy matrices that the special solutions can be written in terms 
of elliptic functions and the normal variational equations are in a type of Lame equations. 



3 Non-integrability theorem 

3.1 Non-integrability proof in the integrable limit 

The chief purpose of the present section is to give a theorem telling the non-integrability 
by using the monodromy matrices obtained in Section 2. We consider monodromy ma- 
trices g. If the eigenvalues {cti, erf 1 , • ■ ■ , cr n , a"" 1 } of monodromy matrices g do not satisfy 
the following relation 

= 1 (3-1) 

for any set of integers {h, ■ •■ ,l n } except the trivial case h — I2 — ■ ■ ■ — l n — 0, we 
call the monodromy matrices g non-resonant. It is already known that the existence of 
a non-resonant monodromy matrix is a basic assumption in order to perform Ziglin's 
analysis JI9"]. Moreover, if there are straight-line solutions such as ( |2 . 1 1| ) whose mon- 
odromy matrices are non-resonant and if the variational equations can be diagonalized 
into a decoupled form like ( |2.33|) , Ziglin's theorem can be generalized to Yoshida's the- 
orem for Hamiltonian systems composed by kinetic energy terms and potential energy 



terms 16 



Yoshida's theorem asserts the following in terms of two different monodromy matrices 
{Oai 9b}' Suppose that there exists an additional complex analytic integral, which is 



holomorphic along the solution (|2.9|), and that one of the monodromy matrices g a is 



non- resonant. Then it is necessary that one of the following two cases, namely, 
(I) <?&(Aj) must preserve the eigendirection of g a {^j) , i-e., g a {^j) must commute with 
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(II) gb{^j) must permute the eigendirection of g a (Xj), i.e., gb{Xj) is written by 
in the base of g a having the eigenvalues i and — i 



p 

-i 



for some suffix j, at least occurs for any other monodromy matrix represented in 
the basis of g a . As for the non-resonance condition of the monodromy matrix g* = 
g\g2g\ g^ , the following lemma is already known: 



Lemma 1 ( ||1 lj , J2[],1994) The n quantities {y25 — 24cos^j\j = 1, • • • ,n} are ratio- 
nally independent. 

We remark that the rational independency of the set 
{^25 — 24cos^j-|j = 1, • • -n} guarantees that the commutator g* = g^g^g^ 1 whose 
eigenvalues are given in Eq. (|2.46| ) is non- resonant. 

For connection with the algebraic number theory on the cyclotomic field Q(exp(^-)) 
over Q see Ref. [|12|]. Thus, we can regard that the non- resonance hypothesis holds for the 
monodromy matrix g*. Now, using Yoshida's theorem and the variational analysis in the 
former section, we obtain a theorem which proves the non-integrability of a FPU lattice 
for an arbitrary set of the system parameters {n, /i 2 , /i 4 |n(> 3) : odd, /i 2 > 0,/i 4 > 0} as 
follows: 



Theorem 1 A FPU lattice ( \2.o) ) which is characterized by an arbitrary set of the system 
parameters {n, fi2, fi^\n > 3,n : odd, \i2 > 0, /14 > 0} has no analytic first integrals besides 
the Hamiltonian itself for fully small energy e(~ 0). 

(Proof of Theorem ELT) 



From the non-resonance hypothesis on g*, we can apply the Yoshida theorem |T6 
to the FPU lattices. According to the above argument, if we prove that at least one 
monodromy matrix g s € {gi,g2} does not have the following properties 

(a) g s (Xj) preserves the eigendirection of g*(\j) 
and 

(b) g s (Xj) permutes the eigendirection of ^(AjXthe eigenvalues of g s {Xj) are i, —i), 

at once for any suffix j(l < n + 1/2) < n) of g s {Xj), then the FPU lattices || are 
concluded to have no other analytic conserved quantities besides the Hamiltonian itself, 
i.e., the assertion of the present theorem holds. In the following, we will show that gi{Xj) 
for any j(^ has neither the property (a) and the property (b). 
When we take the limit e — >• 0, the relations 

74 = /i 4 C 2 , 72m = H2mC 2m - 2 -0, a ^^2, k^0 (3.2) 
hold and we can compute the periods of the monodromy matrices g\, g2 in the limit as 

Ti -> ^==tt, T 2 -> —L=-k + ioo (3.3) 
V 2 72 V 2 72 
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by the formula (|2.27|) . To compute the monodromy matrix gi{Xj) in this limit, we 
rewrite the variational equations ( |2.39|) in terms of the modulus of the elliptic integral k 
as follows 



77T 2k 2 
e ? +4/i 2 sin 2 ( o , \ - )(1 + — C n 2 (k-at))C = 0, j = 1, 



'2(n + l 



1 -2k 2 



n 



(3.4) 



by virtue of equations ( p.!5| ) and fl2.22| ). For arbitrary k G [0, 1], we have a fundamental 
system of solutions as t), £j(k; £)} of (|3.4|) satisfying 



l, |e(^;t)|t=o = e(fc;0) = o 



0) - ^, dtSi , 

g(*;0) = J |g(*;*)|t=o = g(*;0) = i. 



ill 



(3.5) 



Because cn 2 (k; at) in Eqs. 



has the following Taylor expansions JTJ at k = A; 2 = 0: 



cn (k;at) = cos (— t) - -nsin (t^)cos (— t) + 



(3.6) 



and analytic in k at k = 0, the fundamental system of solutions {£ a , £ 6 } are also analytic 
in k at k = || as 



(*l i) = # (t) + $i(t)« + ^ 2 (*)« 2 + 
W*)=W)+d,i(*)* + ^ 2 (*)« a + 

where the unperturbed parts {£| ; o> £j,o} are obtained by 



(3.7) 



^o(t) = cos(2sin(^ iy ) v ^t) 



&>(*) 



-sin(2sin(= ? ^ iy ) > /^t) 



2 v^Sin(^r^-~-^2(n+l) ■ 

from the initial conditions ( |3.5| ). The monodromy matrix gi(Xj) is given by 



(3.8) 



(3.9) 



and with the use of the fact that the fundamental system of solutions (|3.8j ) are analytic 
in k at ft = 0, we show that gi(\j) for j = 1, ■ • • , n have also the Taylor expansions as 

9i(*j) = fl'i.oCAj) + 0i,i(A^)k + ^i,z(Ai)« 2 H • (3.10) 

Thus now, when we consider the low energy limit (e — > 0, k — > 0), we obtain the formula 



-2sin( 



cos 



J 71 " 



(V27rsin(^)) 



// 2 sin(V27rsm( 1 ^ I y)) 



This means that the eigenvalues of gi{Xj) tend to 

jn 



edn( V27rsin( ^^jj) ) 



cos V^TTSin 



{exp (ztt ( v^sin ( — ^— ) ) , exp {-m ( v / 2sin( — 3 -^—- ) ) } 



2(n+l) 



71 



(3.11) 
(3.12) 
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and gi{\j) for any does not have the property of (b). 



Now assume that (^(A?) for some has the property of (a). Then if gi(Xj) as 



well as (72 (Aj) has also the property of (a), we have 

g*(\j) = gi{.X j )g 2 (X j )g^\X j )g 2 \X j ) = id, 



(3.13) 



where id denotes the 2x2 identity matrix. This relation ( [3.13 ) means that gi{Xj) and 
<?2(Aj) commute each other and clear contradicts the non-resonance hypothesis of g*. 
Consider the other case where g\(Xj) has the property of (a) and g 2 (Xj) does not have 
the property of (a). However, in the representation of gi(Xj) and g 2 (Xj) in the basis of 
g*(\j) as 



/i 



i 



) 02 (Aj) 



fl'i(Ay) ; 
the following relation 

9*(Xj) = 9i{X j )g 2 (X j )g^ l (X j )g 2 l (X j ) 



a b 
c d 



(ad — be = 1), 



ad — fi 2 bc ab(fi 2 — 1 
cd(±> - 1) ad bc 



(3.14) 



(3.15) 



must be satisfied. Since g*(Xj) is assumed to have a diagonal representation as g*(Xj) 
-l 



diag 



and from the relation ( 3.15 ), we obtain 

a = 0, d — 0, bc = —1, 



(3.16) 



when g*(\j) 7^ id; g 2 (Xj) must have the property of (b). 
Therefore, in the basis of g*(Xj), we have 



</i(Aj) 



// 
I 



and <?2(Aj 



/3 

-i 



(3.17) 



These relations (|3.17 ) result in 

= g*(Xj) = 9i(X j )g 2 (\ j )gi 1 (\ j )g2 1 (X j ) = g\(X 



*i 
0-7 



.))■ 



(3.18) 



where <jj = —exp{niy25 — 24cos^-}. However, the relation ( pj.!8|) causes again a 
contradiction with the fact that the eigenvalues of gf(Xj) approach to 



{exp(ivr(2v / 2 sin 



J7l 



2(n + l 



■)), exp(— OT(2v / 2sin( 



JTC 



2(n + 1 



"))} 



(3.19) 



in the limit e — > and there is a difference between the eigenvalues of g*{Xj) and the 
eigenvalues of gf{Xj) as will be shown in the Appendix. 

We have seen that gi{Xj) for any has neither the property of (a) or the 

property of (b). Now the theorem holds. (aEnd of proof of Theorem [j]) 
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4 Classification of non-homogeneous nonlinear lat- 
tices 



Here in the present section, we consider the classification of these non-integrable nonlin- 
ear lattices via the degree of non-integrability. We define the degree of non-integrability 
as follows: 

Definition 1 We write a relation 

{H,e}~{H',e'} (4.1) 

if and only if Hamiltonian systems with n degrees of freedom H and H' have the same 
number of additional analytic integrals which are functionally independent together with 
the Hamiltonians H and H' respectively. Here e(e') denotes the total energy of the Hamil- 
tonian systems H(H'). 

If we associate the positive integer 

e = 2 r ■ T~ r ~ l (4.2) 

with a Hamiltonian system H with n degrees of freedom which has r additional analytic 
integrals which are functionally independent together with the Hamiltonian H, we can 
classify Hamiltonian systems via the degree of non-integrability as follows: 

{H, 6} ~ {H\ e'} <=► g({H, e}) = g({H', e'}). (4.3) 

For a non-homogeneous nonlinear lattice 

1 71 771 71 

H^ 4 ,...,, 2 Jq, P ;t) = -EP»+E5E(*-i-*f (4-4) 

z i=i fc=i ZK i=i 

of degree 2m, consider the following transformation 3? from a non-homogeneous nonlinear 
lattice to another non-homogeneous nonlinear lattice: 

Definition 2 We call a transformation 3? from H^ 2tfJj4t ...^ 2m (q, p; t) to H^i ^ (q', p'\ t') 
a homogenizer if the relations 

t' = at, q' = -A-g, p' = -znrP, 

a m—1 a m—L 

a "i-i 

are satisfied where a is a real quantity larger than unity (1 < a < oo). 
Remark 1 

The transformation 9ft in Eqs. Q4.5Q can be regarded as a generalization of the scaling 
transformation for homogeneous nonlinear lattices #j la =o ) -,/*2m-2=o,/i 2m - Moreover, the 
equations of motion itself remain invariant under the transformation 9ft, while a Hamil- 
tonian is changed according to the variations of coupling constants {yU2 7 |l < 7 < to}. 
Thus, we can say that 9ft preserves the degree of non-integrability as 

{9WJ,9fc} ~ {H,e}. (4.6) 



'm—1 



/ 1 2m 



l^2r, 



(4.5) 
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Remark 2 

The operation 3ft corresponds to a kind of renormalization group along the time axis 
from the view point of the coarse graining in statistical physics and it is easy to see that 
the fixed point of 3ft is the homogeneous nonlinear lattice H*, i.e., 

h* = mr, (4.7) 

where 

_ 1 v^n 2 i M2m v^n I „ _ _ \2m 
— 2 l^i=\Pi 2m 2^i=x\Hi-l Hi) 

and 3ft' denotes Z iterations of the operation 3ft. This is the reason why the 3ft is called a 
homogenizer here. 

Remark 3 

Under the transformation 3ft, the total energy e is changed into e' by the formula 

e' = 3fte = -ire. (4.9) 



By the above remarks, we can conclude that 



{H^ 2 ^-^e}~{H^ 4 ,e} (4.10) 



for an arbitrary FPU lattice {if M2iAt4 , e}. If we take the limit a 1 — ► oo preserving the re- 
lation -Ij- « 0(1), L.H.S. of fl4.10p approaches a homogeneous nonlinear lattice of degree 
4, which means that the degree of non-integrability of a homogeneous nonlinear lattice 
is the same as the degree of non-integrability of a FPU lattice in the high energy limit 
e — > oo. In relation to the non-integrability of the homogeneous nonlinear lattice, there 
is a known result by Yoshida|]l7| stating as follows: 



Considering n number of distinct exponents 
A 



x 3 



8m(2m — 1) Q . 77T , . 

+ ^ ^sin 2 — ^— - for j = l,.-.,n), 4.11 

or the n representative of the Kowalevski exponents^L^ of the homogeneous nonlinear 



lattice of degree 2m. If the n quantities {A^} in ( 4.11 ) are rationally independent, then 
the homogeneous nonlinear lattices have no analytic first integrals except the Hamilto- 
nian itself for 2m(> 4) and for odd n(> 3). 



Note that the n quantities ( 4.11] ) equal the phase factors ( |2.49| ) of the eigenvalues of 



the commutator g*. Thus, the non-resonance hypothesis proving the non-integrability in 
the low energy limit e — > induces also the non-integrability of nonlinear lattices with 
a nonhomogeneous potential function in the high energy limit e — ► oo. The reason for 
this uniqueness of the non-resonance hypothesis lies in the fact that after a conformal 
mapping z = 4 (t), the closed loop of the commutator g\g29\ X (h} becomes a homotopy 
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equivalent loop which corresponds to Yoshida's construction of the monodromy matrix 
G — (GqG% ) 16 giving the same non-resonance hypothesis in terms of Kowalevski expo- 
nents. Here, G corresponds to a counterclockwise closed loop around the singularity 
z = and G\ corresponds to another counterclockwise closed loop around the singular- 
ity z = 1 in the complex z plane, where z satisfies the following Gauss hyper geometric 
equation: 



z(l 



dz 2 



+ 



d? 3 , 
Tz + 2 74Sm 



JTT 



-X = o 



(4.12) 



'2(71+1)' 

into which the conformal mapping z = <f) A transforms the variational Eq. ( |2.33 ) with the 
condition 72 = and Go, G\ are the two fundamental monodromy groups of the Gauss 
hyper geometric equation. See also Fig. 3. 



So far, we have succeeded in checking the non-integrability of the FPU lattices both in 
the low energy limit and in the high energy limit from a single non-resonance hypothesis. 
However, there still remain some problems about the integrability of FPU lattices for a 
finite energy < e < 00. To consider the problem, we introduce a real parameter 



X 



€fJL4 



> 0, 



(n + 1)m| 

where x is found to be a dimensionless parameter, because 



(4.13) 



Dimension of x 



\ L 2 ' 




1 


Lmt". 




Lmtt/J 


[1] 






_M 2 T 4 _ 





[1] 



(4.14) 



by the standard unit system as [M] = [kg] , [T] = [sec] , [L] = [m]. We remark here that 
X is invariant under the transformation 3? preserving the degree of non-integrability. 
This means that it is sufficient for us to classify the degree of non-integrability of the 
FPU lattices by using this dimensionless parameter x- Furthermore, with the use of the 
parameter x, we can replace the formula of the modulus of the elliptic function k in Eq. 
(PI) by 

' ' - ' (4-15) 



v/2V VT+Ix' 

We note here that k — > if and only if x ~~ *■ and this limit can be realized for each set 
of {e, /X2, A*4|0 < e, /X2, A*4 < 00} when we take n(:odd) a fully great number. If we recall 
that in Theorem 1 it is essential for proving the non-integrability of the FPU lattices to 
take the limit k — > 0(x — > 0) so that the eigenvalues of gi(Xj) approach the quantities in 
Eq. Q3.12 ), we have the following theorem: 



Theorem 2 A FPU lattice ( \27q ) characterized by an arbitrary set of the system parame- 
ters {e, /X2,A t 4|0 < e < oo,0 < fi2 < oo,0 < fi4 < 00} has no additional analytic integrals 
of motion besides the Hamiltonian itself for a fully great number of degrees of freedom 
n(n : odd). 
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This result shows that we need a fully great number n of degrees of freedom to 
guarantee the strongest non-integrability g({H FPU , e}) = 3 n_1 of each FPU lattice H FPU 
with the systems parameters {/i 2 ,/-t4|0 < /i 2 < oo, < /i 4 < 00} for < e < 00 and a 
problem remains still open as to whether we can prove the non-integrability of each FPU 
lattices with the system parameters {/z 2 , A*4, n|0 < /i 2 < 00, < fi^ < 00, n : odd, n > 3} 
for a given finite energy e(0 < e < 00). Clearly, Theorem 1 and Theorem 2 have the 
reciprocal conditions in terms of the system parameters {e, /z 2 , /Z4, n} and we can unify 
them into a single theorem as follows: 



Theorem 3 A FPU lattice characterized by the system parameters 

{e, /x 2 , /i 4 , n|0 < e < 00, < yU 2 < 00, < yU 4 < 00, n > 3, n : odd} has no additional 

analytic integrals of motion besides the Hamiltonian itself when the following condition 



(n + !)//! 



zs satisfied. 



5 Summary and discussion 

We have got the non-integrability proof of the FPU lattices by tracing the following steps. 
In Section 2, we have shown that the singularity analysis can be performed towards the 
non-homogeneous nonlinear lattices by finding the special straight-line solutions in terms 
of elliptic functions. Especially, it is shown there that the normal variational equations 
can be decoupled into separated equations and each of them happens to be a Lame 
equation and consequently we can compute the eigenvalues of the monodromy matrices 
associated with the counterclockwise loop along the period cell of the special solutions. 
In Section 3, with the use of Yoshida's theorem whose validity is guaranteed when that 
the normal variational equations can be decoupled, we can prove the non-existence of the 
additional analytic conserved quantities besides the Hamiltonian itself for FPU lattices 
in the low energy limit. In Section 4, we have considered the classification of FPU lattices 
via the degree of non-integrability. There, by introducing a transformation 3? from a non- 
homogeneous nonlinear lattice into another non-homogeneous nonlinear lattice, which 
preserves the equations of motion, we have shown that the degree of non-integrability 
of each FPU lattice in the high energy limit is the same as that of the homogeneous 
nonlinear lattice whose non-integrability is known to be proven by the analysis using 
Kowalevski exponents. In case of FPU lattices with arbitrary energy in < e < 00, 
we have a theorem showing the non-integrability for a fully great number of degrees 
of freedom by introducing a dimensionless parameter which characterizes the normal 
variational equations. This theorem (Theorem 2) and Theorem 1 have the reciprocal 
sufficient conditions for the non-integrability of FPU lattices and finally they are unified 
into a single theorem (Theorem 3) by the dimensionless parameter \. 

There remains an open problem to be solved: Can we have a theorem about the 
non-integrability of more general nonlinear lattices of degree 2m > 4? More specifically 
speaking, can we have a theorem asserting the non-integrability of the non-homogeneous 
nonlinear lattices of degree 2m > 4 such as Theorem 1 and Theorem 2? The present 
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neck is in a fact that the normal variational equations along the special solutions of 
hyper- elliptic functions do not belong to a class of Lame equations, though they can 
always be decoupled into n — 1 separated variational equations, and we cannot compute 
the eigenvalues of the monodromy matrices associated with these variational equations 
in general and therefore we cannot perform the Ziglin analysis completely. However, the 
problem on the degree of non-integrability for more general non-homogeneous nonlinear 
lattices remains an interesting open question to be solved, because the present analysis 
on the non-integrability of the FPU lattices, together with Yoshida's argument about 
the homogeneous nonlinear lattice in terms of Kowalevski exponents, strongly suggests 
universality about the non-integrability of more general non-homogeneous nonlinear lat- 
tices in Eq. (4.4). 



6 Appendix 

In this appendix, we show the difference between the eigenvalues of g*{\j) and the 
eigenvalues of gf(\j) for any Sf ^)- We prove this by reductio ad absurdum. First, 
we assume that the following equality 



SpecG^A,)) = {exp(z7r(2^sin(^)),exp(-m(2y^sin( 1I ^ Ty ))} 



{-exp(7n J25 - 24cos^ T ), -exp(-W25 - 24cos^)} = Spec(^(A,)) 



(6.1) 



would hold for some Our purpose here is to show that the relation ( |6.1|) is 

false. By (|6.1|), one of the following equalities 



'25 - 24cos^- + 1 + 2m+ = 2 v / 2sin( - ^ - ) (6.2) 
n+1 + K 2(n + l) J V ; 



J 25 - 24cos-^- + 1 + 2m_ = -2 v / 2sin( , ^ A (6.3) 
V n + 1 y 2(n+l)' v ' 

holds for m + , m_ G Z and j ^ Because of the inequalities 



1 < W25 - 24cos^- < 7, < 2 v / 2sin(— ] -^— ) < 2y/2, (6.4) 

y 71/ ~\~ ±. 2ii Tt ~~\~ 1 ) 

it is sufficient to consider the equalities (|6.2|) and ( |6.3| ) for m + e {—3,-2,-1,0} and 

m_ G {—5, —4, —3, —2}. If we rewrite sm( 2 (n+i) ) as X, Eq. (|6.2j ) and Eq. ( |6.3|) are 
changed into the following equations 

AQXl + (1 + 2m + )4 v / 2X + + 1 - (1 + 2m + f = 0, (6.5) 

40X 2 - (1 + 2m_)4 v / 2X_ + 1 - (1 + 2m_) 2 = (6.6) 
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respectively. By considering X( 
solutions as follows: 



X + , X_ > 0, we can easily obtain a series of the 



X+ 

x + 

x + 
x_ 

x_ 



(m + 
(m + 

(m + 
(m_ 
(to_ 
(to_ 



5^4 



0) = 

-1) 

-2) 
-3) 

6 ) - 10 



20 



ggO > 7+V289 



20 



(6.7) 
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Furthermore, because of the following relations < X± < 1, X± ^ 4g for j 7^ we 
can discard the solutions X + {m + = 0),X + (m + = — 2),X + (m + = — 3),X_(m_ = —4). 
For the solutions X + (m + = — l),X_(m_ = —2), we have the equalities as 



cos 
cos 



J7T 

n+1 
n+1 



l-2(X + (m + = -l)) 2 
1 -2(X_(m_ = -2)) 2 



24 
25 
21 
25 



respectively. However, these relations 



cos- 



G Q 



cos- 



1|) contradict with a fact [PJH that 

J7T 



n+1 n+1 

Thus, we can also discard the solutions X + \ 
consider the solutions X- (m_ = 
following relation 

25C 4 

is satisfied for the solution X_ ( 



0. 



or 



1 
2' 



or 



1 
2" 



(6.8) 



(6.9) 



-3) and X. 



— 1) and X_(m_ = —2). Let us 
-5). If we rewrite exp^y as £, the 



70C 3 - 46C 2 + 70C + 25 = 
,_ = —3), while the relation 



25C + 462C 3 + 626C 2 + 462£ + 25 = 



(6.10) 



(6.11) 



is satisfied for the solution X_(m_ = —5). Since ( is one of 2(n + l)-th root of unity 
and 2(n + 1) (mod 4) = 0, the minimal polynomial -P(C) with coefficients of rational 
integers which has a solution ( for P(() = must be one of the cyclotomic polynomials 
$i(y) G Z[Y\: 

<p(i) 

HY) = l[( Y -0), (6.12) 

i=i 

where 2(n + 1)( mod /) = 0, {d} is a set of all primitive I th roots of unity and ip{l) is 
the number of positive integers which are less than or equal to / and relatively prime to 
/ and this function tp is called Euler's ^-function. All cyclotomic polynomials <&i(Y) are 
known to be irreducible over Q and we can check that cyclotomic polynomials $i(Y) 
whose degrees are less than or equal to 4 are restricted to the case[T^] that 



N = {&i(Y)\2 < Deg$i(y) < 4} 

= {$ 3 (r) = Y 2 + Y + 1, $t{Y) = Y 2 + 1, $ 5 (F) = F 4 + Y 3 + Y 2 + Y + 1, 

$ 6 (V) = y 2 - y + 1, $ 8 (y) = r 4 + 1, $ 12 (y) = r 4 - r 2 + 1}. 



(6.13) 
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Using this fact, we can say that a cyclotomic polynomial $i(Y) £ K must divide either 
the polynomial 

P X {Y) = 25Y 4 + 70Y 3 - AQY 2 + 70Y + 25 (6.14) 
in Eq.( |6.10| ) or the polynomial 

P 2 (Y) = 25F 4 + 462F 3 + 626F 2 + 462Y + 25 (6.15) 

in Eq. fl6.lip . It is easy to check that this is impossible. Thus now, we can conclude that 
the assumption fl6.1|) is false; i.e., we have shown that for any j(^ 

Spec(^(A,))^Spec(^(A J )). (6.16) 
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